A stellar test of the physics of unification 
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We discuss the feasibility of using solar-type main-sequence stars as probes of fundamental physics 
and unification. We use a simple polytropic stellar structure model and study its sensitivity to 
variations of the gravitational, strong and electroweak coupling constants in the context of unification 
scenarios. We quantify the sensitivity of the Sun's interior temperature to these variations, finding 
|Aa/o| < 1.3 x 10 -4 for a 'canonical' choice of unification scenario, and discuss prospects for future 
improvements. 
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I. INTRODUCTION 

Cosmology and particle physics are presently experi- 
encing a truly exciting period. On the one hand, both 
have remarkably successful standard models which 
are in agreement with a plethora of experimental and 
observational data. On the other hand, there are also 
strong hints that neither of these models is complete, 
and that new physics may be there, within the reach of 
the next generation of probes. 

In this context, it is important to identify laboratory or 
astrophysical measurements that can give us more infor- 
mation about the nature and properties of this still un- 
known physics. Here we will discuss one such example. 
Several authors have suggested that stellar physics can 
be used to test the behavior of the gravitational sector 
@~Bl- We will go further, as in by showing how stel- 
lar physics can also provide information on the physics 
of unification and on the possible spacetime variation of 
nature's fundamental couplings. 

Nature is characterized by a number of physical laws 
and fundamental dimensionless couplings, which histori- 
cally we have assumed to be spacetime-invariant. For the 
former this assumption is a cornerstone of the scientific 
method, but for the latter it is an assumption with no 
justification @, 0]- There is ample experimental evidence 
showing that fundamental couplings run with energy, and 
many particle physics and cosmology models suggest that 
they should also roll with time [8j. This happens in cos- 
mological models with dynamical scalar fields, including 
string theory. 

Do the fundamental constants vary? Answering this 
question has key implications for cosmology and funda- 
mental physics, and in particular can shed light on the 
enigma of dark energy [j| [lOj . The current observational 
status is controversial, with claims for and against 



variations of the fine-structure constant a and the 
proton-to-electron mass ratio n at redshifts z ^ 1 — 3. As 
shown by recent workshops devoted to the topic the de- 
bate is open 17J and a resolution demands better data, 
but also independent ways to search for these variations, 
which may confirm or contradict these indications. 

In addition to a and /x, the gravitational constant G 
is also interesting to study. Speaking of variations of 
dimensional constants has no physical significance: one 
can concoct any variation by defining appropriate units 
of length, time and energy. However, one can always 
choose an arbitrary dimensionfull unit as a standard and 
compare it with other quantities. If one assumes particle 
masses to be constant then constraints on G are in fact 
constraining the (dimensionless) product of G and the 
nucleon mass squared. A better route is to compare the 
strong interaction with the gravitational one: this can 
be done by assuming a fixed energy scale for Quantum 
Chromodynamics (QCD) and allowing a varying G, or 
vice-versa. With these caveats, probes of a rolling G 
provide key information on the gravitational sector [l8j | - 

Paradoxically, G was one of the first constants to be 
measured (by Cavendish) but is now the least well known, 
a result of the weakness of gravity. In fact it is not known 
to better than percent level |.19]. While a rolling a or /i 
implies a violation of the Einstein Equivalence Principle, 
a rolling G is compatible with it, but it does violate the 
Strong Equivalence Principle. Any Grand-Unified The- 
ory predicts a specific relation between variations of a 
and \x [2fj| , and simultaneous measurements of both pro- 
vide key consistency tests 21] . Different classes of grav- 
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itational theories also predict different relations between 
these variations and those of G, so tighter gravitational 
constraints can discriminate between them. 

Our work, which is in the same spirit as [!, [23|, 
considers the feasibility of using stars as probes of fun- 
damental physics and unification. Specifically, we are 
interested in the impact of possible spacetime variations 
of a, fi and G on the structure and evolution of stars. 
Stars are interesting for this purpose both because they 
are in principle sensitive to variations of all three relevant 
couplings and because their long lifetimes may make even 
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tiny variations noticeable. 

Presently we will limit ourselves to a study of poly- 
tropic stars. Although they are of limited use for de- 
scribing realistic stars, they offer the advantage of math- 
ematical simplicity, which in turn allows one to develop a 
clear physical picture of how the changes in the parame- 
ters impact the structure of the star. In subsequent work 
we will discuss more realistic stars, whose study requires 
detailed stellar evolution codes. 



Similarly for the mass difference between neutrons and 
protons, a = m n ~ m p , we hnd 



^ = [0.1 + 0.7 5-0.6 R}^. 
a a 

For the ratio r\ = m n /m pi we have 



(6) 



^=(i-lW-0.5S+1.4JJ3^. (7) 



II. PHENOMENOLOGY OF UNIFICATION 

We wish to describe phenomenologically a broad class 
of models which allow for simultaneous variations of sev- 
eral fundamental couplings, such as the fine-structure 
constant a = e 2 /(he), the proton-to-electron mass ratio 
fi = m p /m e and Newton's gravitational constant G. The 
simplest way to do this is to relate the various changes 
to those of a particular dimensionless coupling, typically 
a. Then, if a = ao(l + S a ) and 



AA _ fc Aq 
A a 



(1) 



we have A = Aq(1 + k,A&A) and so forth. 

The relations between the couplings are model- 
dependent. Here we will follow [20j], considering a class of 
grand unification models in which the weak scale is deter- 
mined by dimensional transmutation and further assum- 
ing that relative variation of all the Yukawa couplings is 
the same. Finally we assume that the variation of the 
couplings is driven by a dilaton-type scalar field (2~3 |. In 
this case one find that the variations of p and a are re- 
lated through 



^ = [0.8 R - 0.3 (1 + 5)] — 
fx a 



(2) 



where R and 5 can be taken as free phenomenological 
parameters. Their absolute value can be anything from 
order unity to several hundreds, but while R can be posi- 
tive or negative (with the former case being more likely), 
physically one expects that 5 > 0. Nevertheless, we can 
simply treat both as phenomenological parameters to be 
constrained by astrophysical data. 

Concerning Newton's constant G, it is easiest to work 
in terms of the dimensionless couplings 
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Gmj 
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(3) 



For our purposes it's natural to assume that the QCD 
scale and particle masses vary, while the Planck mass is 
fixed. We then have 



Am. 



(1 + 5) 



Aa 



and 



2^ = 2 [0.8i? + 0.2(l + 5)l — 

m„ a 



(4) 



(5) 



Relative variations of quantities such as the neutron life- 
time and the deuteron binding energy can also be cast in 
this form 20]. 

One important effect of varying masses is that they 
lead to a change in the star's density, even though the 
number density has the standard behavior. Indeed 



Am 



Aa 



Pstd 



1 + ^ = 1 + [0.8 # + 0.2(1 + 5)] — ; (8) 



thus the importance of this effect depends not only on 
the the value of the a variation but also on the values of 
the phenomenological parameters R and 5. 

If one takes at face value the current evidence for a 
varying a [HI [HJ and the null results for p [ll| [ll| , one 
could infer that 



Au Aa 

— < 

fi a 



(9) 



If so, we can use the above fi-a drift equation to obtain 
the following (approximate) relation 



#~-(l + 5) 



(10) 



which can in principle be used to simplify the analysis: 
in this case we have only two additional free parame- 
ters, Aa/a and R. Although our analysis below will be 
generic, we will comment further on this particular limit. 
There are also [13j indications of a spatial a dipole, in 
the direction RA = (17.3 ± 1.0)h, S = (-61 ± 10)°. If 
confirmed, this would suggest that stars in different di- 
rections would be subject to different couplings. 



III. POLYTROPIC STARS 

A polytropic star (see [25[) is a simplified model for the 
structure of a star in equilibrium, built from the equation 
of mass and assuming hydrostatic equilibrium. Specifi- 
cally, it uses 



dm 
dr 



Airr 2 p 



and 



dp 
dr 



Gmp 



(11) 



Here r is the radial distance to the centre of the star, 
m the mass within the sphere of radius r, p the density, 
and p the pressure. These equations can be solved if 
an equation of state of the type p = p(p) is set. The 
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polytropic solution corresponds to a symplified equation 
of state in the form 



P 



= K p 1 



-l/r. 



(12) 



where K is the polytropic constant (to be defined by the 
boundary conditions) and n is the polytropic index (to be 
selected). The use of this relation avoids the need to in- 
clude an additional equation for temperature, describing 
the transport of energy in the interior. 

In such a simplified approach the structure of the star 
is described by the numerical solution of the Lane-Emden 
Equation, as given by 



_f_ _d / 2 dfl 



(13) 



In order to obtain this equation from Eqs (|TTj) . the fol- 
lowing variables have been used; 9 n = pj p c and £ = r/a, 
with p c being the value of the density at the centre of the 
star and 



K(n+1) 
AnGp 



l-l/n 



(14) 



The initial conditions, at the centre, for the Lane-Emden 
equation are 9{t\ = 0) = 1 and 0'(£ = 0) = 0. The 
solution goes from £ = to £ = £ s , where 8(£ s ) = and 
9' s = g|(£ s ) < 0. This value of £ s , corresponds to the 
first zero of the solution for 9 for each value of n used. 
There is a finite value of £ s > y/6 for any < n < 5. 

If we consider a star with total mass M* and radius i?* 
(not to be confused with the unification parameter R), 
then 



Pc = 



39' j AirRl ' 



and 



K 



G 
n+1 



4tt 



MT 



-il/n 



n+1/ 



(15) 



(16) 



The production of energy is not required, as any solu- 
tion can be found with the selected radius and mass for 
a value of the index n £ [0, 5[. However, if the luminosity 
of the star (£*) is known, it may be used as a boundary 



condition that can establish the value of n. Finally, if this 
boundary condition for the energy is used, then we can 
obtain the polytropic index n. This is done by imposing 
that the luminosity of the polytrope is calculated from 



M 



(17) 



where e is the emissivity, or energy production rate per 
unit mass and time, inside the star. The expression for 
e will depend on the process used by the star to produce 
energy in the interior. Therefore it will be given by the 
physics describing the fusion of hydrogen into helium in 
main sequence stars. The version we adopt here are the 
global fits for the emissivity of the PP Chains (relevant 
for low mass stars, like the Sun or late- type stars). 

When doing so, we need to assume an equation that 
provides the values of the temperature from density and 
pressure, for a known chemical composition. For the 
present work the ideal gas law is used, being an ade- 
quate approximation for stars, similar to the Sun, in the 
main sequence. The value of L+ provides the condition 
to select n, if we require that the condition L n = L* is 
satisfied. 



IV. IMPACT ON THE STAR FROM a 

Since the variation of fundamental couplings directly 
affects the units we use to measure quantities with dimen- 
sions, it is convenient to work with dimensionless quan- 
tities. We will therefore use the dimensionless version of 
the previous equations by using the relations; 



— , j— — -, p^f 
m p trip n L c~ L 



Thus for density we have 

Pc _ M+frrip 



P 



3„3 



rripH J ' 



4tt(91 



(18) 



(19) 



while for the polytropic constant 



K 



m p i/n frV™ c 2 " 3 /" 



(AU/m p ) n 
_{R i ,/rrip 1 K L c- 1 



n — 3 



1/n 



47T 

a +1 (-8' s y 



l/n 



(20) 



For the emissivity, £y, of the reaction between species i and j, the following family of expressions can be used (see 
Adelberger et al. [26[ for details, and references therein); 



3Ai 



2 \ 1/6 



So K-^p 1 - 2 ^ 71 a 1 / 3 



m^h 2 m p h- 3 c 5 / 3 6»"+ 2 /3 



exp 



1/3 



a 



2/3 



T ^ l/n n 

Rpc c- 2 



31/3 



(21) 
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where / is the energetic gain of the nuclear reaction in units of m p c 2 , Aij is the reduced atomic mass in units of m p , 
Wij = ZiZjAij where z is the atomic number and Sq the reference value for the S-factor (the effective cross section). 
By using an expression for the emissivity, based on the dominant reaction of the PP Chain, we may then find the 
relation based on luminosity, that provides the relation that determines the value of n. This is, 

_ 3 



(k p s *y 

-2 fc 3 3 

TOr, ft2C2 



fa 



a i, 



i2an 



mlh^c 4 



(22) 



A model of a specific star, when we have (-M*, i?*, L*), 
can be obtained by solving Eq. (TlUl) together with 
Eq. (j2"2"|) , where we impose that L n = L*. The solu- 
tion provides the behavior of p and p. Any change in the 
physics will imply a different solution in the interior for 
the same global parameters. The model accommodates 
such a change by having a slightly different value of n 
that fulfills the equation for the luminosity. 

This approach allows us to obtain quantitative esti- 
mates of the relative change of the stellar structure when 
the underlying physics changes. Although a polytrope is 
not the best available model of a real star, it is a close 
enough model for our purposes. The assumption behind 
this approach is that the relative change from the real 
star is of the same magnitude as the relative change of 
the polytropic model, when every other unknown aspect 
of the physics is not allowed to change. 



V. RESULTS AND CONSTRAINTS 
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FIG. 1: The relative change of the Sun's central temperature, 
for Aa/a = — 10 -4 , as a function of the unification parameter 
R (with the choice S = 0). Allowing a maximal variation of 
one percent yields the limit \R\ < 90. 



From observations of a specific star, the value of its 
stellar parameters may be obtained, specifically its mass 
ML , radius i?* and luminosity L+ (e.g. Creevey et al. 
[27j |. and references therein). This allows us to build a 
model, assuming that the set of equations and physics 
we use are valid, that fits these boundary conditions and 
describe the structure of the interior. However if our 
physics is changed by a small amount, we must find a 
new model, adjusted to the same boundary conditions, 
but with a slightly different structure. This approach can 
also be implemented for polytropic models, which allow 
us to quantify with great precision the change that one 
single aspect of the physics can have on the structure 
of the interior. Specifically, we can determine by how 
much specific thermodynamic quantities will change in 
the interior when we introduce a small change in the 
physical constants entering the basic equations of stellar 
structure. 

In order to have a useful test, we must refer to a quan- 
tity that can be obtained from the observation of real 
stars. The ideal candidate is the value of the tempera- 
ture near the centre, since this value has a direct effect on 
sound speed or neutrino production. The former can be 
expected to be measured with significant precision with 
asteroseismology (or helioseismology for the solar case) 



while the latter may be used for the Sun. We do not in- 
clude the accumulated effect in time, since that requires 
a proper modeling of the evolution which is beyond the 
scope of the present work. 

Presently, the precision we have for the temperature 
(or sound speed) in the interior of the Sun, is well below 
1% (see [23]). While far from this, a very high precision 
is also expected from the asteroseismology of other stars, 
as is currently being obtained with data from the NASA 
space mission Kepler (e.g. Chaplin et al. [29]). 

We thus use the formalism in the previous sections 
to estimate the impact that variations of a will have 
in the internal structure of a star for a given unifica- 
tion model (parametrized by specific values of R and 
S). We consider a polytropic star similar to our Sun 
(assumed to have a politropic index no), to which will 
correspond a slightly different value for n. This star 
will have the same value of M*/m p , R*/ '(m~ foe \ and 
_L*/(mpS _1 c ) (M*, i?* and L* are the Sun's mass, radius 
and luminosity) but for a range of values where a differs 
from the standard value. We then calculate by how much 
its central temperature (in fact the ratio KbT c / (m p c 2 ), 
where T c is the central temperature) must differ from 
that of our Sun, by solving Eq. (|2"21 to determine the 
value of n. 
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suggested in [20] of R 
\Aa/a\ < 1.3 x 1CT 4 . 



30 and S ~ 160, we find 



FIG. 2: Bounds on the (a, R) parameter space, for the case 
5 = 0. The region above and to the right of the dashed 
line is excluded by our analysis. The dots (crosses) show the 
modulus of the estimated values of the maximum (minimum) 
allowed R for particular choices of a. 



We have fixed specific values of Aa/a and S = and 
let R vary freely. Since we are only considering the PP 
chains the main sensitivity of the polytropic model will 
be to changes in a p (as well as a itself), and therefore 
we can only hope to constrain a linear combination of R 
and S. We may then assume S = and get a constraint 
for R, which from Eq. ((SJ) will in fact be a constraint on 
R + S/A. We have assumed R m in — —Rmax', in the range 
we have considered, 10~ 8 < |Aa/a| < 10~ 2 , this is an 
excellent approximation except at the top limit (where it 
is still adequate although no longer plausible). 

In Fig. Q]we show an example for Aa/a = — 10~ 4 , in 
which case R m ax — 90. Repeating this analysis for differ- 
ent values of a we can identify the region in the {R, S, a) 
parameter-space which is consistent with the currently 
estimated experimental uncertainty for the central tem- 
perature of our Sun. This is summarized in Fig. [3J and 
yields our final bound 



|4R + S| < 1(T 



Ac 



(23) 



As a simple illustration, if we take the typical values 



VI. CONCLUSIONS 

We have used a simple polytropic model of stellar 
structure to quantify the impact of variations of the grav- 
itational, strong and electroweak coupling constants, in 
the context of unification scenarios, on solar- type stars. 
Although our analysis involves a number of simplifying 
assumptions, it serves to show that stars can yield tight 
constraints on these scenarios, with the added advantage 
that they nicely complement other astrophysical probes 
of these variations (such as quasar absorption systems 
and the cosmic microwave background) in terms of the 
characteristic density, length and time scales that they 
are sensitive to. 

Naturally a polytropic model is of limited use for de- 
scribing realistic stars. Nevertheless, of particular rele- 
vance are stars near the lower limit for stellar mass, cor- 
responding to fully convective stars (closely represented 
by a polytropic model with n = 1.5) and having very 
long lifetimes in the main-sequence. Asteroseismology of 
such stars may be able to provide robust limits on the 
central temperature, with this case corresponding very 
closely to the underlying assumptions used in this work. 

The logical next step is to study more realistic models 
of stars, for which one requires detailed stellar evolution 
codes. The very long lifetime of a star (particularly low- 
mass ones) is expected to be a key advantage since even 
tiny effects, accumulated over a long period, can have 
distinct observational consequences. 
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